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LJT  ,  _)  t  ABSTRACT 


An  approximate  analytical  queuing  model  for  the  intraship  transfer 
of  material  phase  of  the  replenishment  at  sea  problem  is  formulated. 
The  physical  model  has  two  parallel  sets  of  three  service  facilities 
in  series  that  move  material  to  a  final  stage  of  three  service 
facilities  in  parallel.   Exponential  service  times  and  finite  queue 
space  are  assumed  for  each  facility,  except  that  the  initial  facilities 
always  have  infinite  queues.   This  implies  that  the  output  from  an 
individual  facility  is  approximately  Poisson.   The  mean  output  rate  for 
the  queuing  system  is  obtained.   Several  parametric  studies  on  critical 
parameters  are  performed,  with  the  results  presented  in  graphical  format, 
The  model  is  easily  enlarged  to  allow  a  general  number  of  facilities  in 
series  in  any  of  the  parallel  paths.   Generalization  to  more  than  two 
paths  in  parallel  for  the  initial  stages  may  be  possible. 
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I.   INTRODUCTION 

Underway  replenishment  (UNREP)  of  the  ships  of  the  United  States 
Navy  at  sea  is  a  mandatory  requirement  in  maintaining  a  continuous 
posture  of  readiness  to  meet  the  present  and  future  global  commitments 
and  responsibilities  of  this  country.   However,  replenishment  of  ships 
at  sea  should  detract  as  little  as  possible  from  the  primary  missions 
of  a  fleet.   This  implies  that  the  goal  of  the  UNREP  operation  is 
complete  resupply  of  combatant  ships  in  the  minimum  possible  time. 

There  are  several  phases  of  the  replenishment  operation,   possibly 
not  all  distinct.   One  list  of  phases  might  be:   Long  range  planning 
required  to  schedule  a  replenishment  both  in  terms  of  time  and  place; 
the  actual  movement  of  ships  toward  the  prescheduled  rendezvous;  when 
several  supply  and  combatant  ships  are  involved  in  a  replenishment 
operation,  the  scheduling  and  movement  of  all  the  combatant  ships 
through  the  replenishment  cycle;  and  the  actual  transfer  of  material 
from  the  hold  of  the  supply  ship  to  designated  storerooms  on  the 
combatant.   It  is  obvious  that  all  phases  require  careful  attention 
if  the  operation  is  to  be  efficient.   Consequently,  all  phases  are 
fruitful  study  areas  in  a  continuing  effort  to  improve  UNREP  efficiency, 

The  study  reported  herein  dealt  with  the  actual  transfer  operation. 
Perhaps  it  is  the  most  important  phase  in  that  the  ships  involved  are 
most  vulnerable,  their  capability  to  maneuver  reduced  and  their 


Throughout  the  paper  the  terms  UNREP  and  replenishment  operation 
will  imply  the  entire  operation;  planning,  scheduling,  shipmovement, 
etc.,  while  the  term  transfer  operation  denotes  the  movement  of  material 
from  hold  of  the  supply  ship  to  storeroom  of  the  combatant. 


watertight  boundaries  open  to  allow  passage  of  material.   Unfortunately 
in  the  past  it  appears  that  this  phase  has  received  less  attention  than 
it  warrants  both  from  the  standpoint  of  detailed  study  and  the  introduc- 
tion of  improved  material  handling  equipment  and  techniques. 

Recent  efforts  have  been  made  that  should  considerably  improve  the 
overall  speed  and  efficiency  of  the  transfer  operation.   Improvements 
include  the  introduction  of  new  supply  ships  with  the  capability  to 
employ  the  unit  load  system  of  material  storage  and  transfer,  more  and 
better  mechanized  material  handling  equipment  on  both  supply  and 
combatant  ships,  and  a  realization  of  the  problem  by  ship  designers 
which  has  resulted  in  improved  match  of  replenishment  stations  between 
supply  and  combatant  ships.   The  unit  load  concept  provides  for  mate- 
rial to  be  stored  on  the  supply  ship  in  the  configuration  for  future 
transfer.   Full  implementation  of  the  concept  will  maintain  the  unit 
load  intact  through  ultimate  storage  in  the  storeroom  or  magazine  of 
the  combatant  receiver.   Although  this  concept  may  never  be  completely 
feasible,  especially  in  the  case  of  smaller  combatant  ships,  it  should 
provide  a  vast  improvement  to  perhaps  the  most  important  and  demanding 
of  all  UNREP  operation^,  the  rearming  of  aircraft  carrier 

The  transfer  operation  can  be  represented  as  a  queuing  problem. 
It  is  a  system  of  parallel  queues  in  series  with  movement  between  the 
parallel  paths.   Attempts  at  reaching  an  analytic  solution  are  further 
thwarted  by  the  finite  waiting  room  in  the  queues  involved.   Formal 
efforts  to  analyze  the  problem  have  been  few  in  number.   The  powerful 
computer  simulation  tool  has  been  the  general  mode  of  analysis.   An 
early  simulation  was  conducted  by  Stanford  Research  Institute  [s] .      A 
detailed  computer  simulation  study  has  also  been  done  by  Alfred  M. 
Feiler  of  Project  TRANSIM,  University  of  California  at  Los  Angeles, 
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Department  of  Engineering,,   Computer  simulation  studies  are  currently 
being  done  at  Stanford  Research  Institute  and  Center  for  Naval 
Analyses . 

II.   BACKGROUND 

A,   THE  TRANSFER  OPERATION 

Although  some  transfer  of  material  between  ships  at  sea  is 
presently  accomplished  by  helicopter,  the  vast  majority  of  material 
moves  between  ships  via  attached  horizontal  rigs.   A  typical  move- 
ment sequence  for  material  transferred  at  sea  would  commence  with  a 
vertical  lift  from  the  hold  of  the  supply  ship  to  the  main  or  transfer 
deck  by  means  of  an  elevator  or  a  boom  and  windlass.   Once  on  the 
transfer  deck  the  material  will  be  moved  to  the  transfer  station  by 
fork  lift  truck  or  hand  pallet  truck.   The  intraship  transfer  of 
material  is  accomplished  by  one  of  several  different  cable  and  pulley 
arrangements.   The  type  employed  on  a  given  transfer  will  depend  upon 
the  types  of  ships  involved  and  the  amount  and  type  of  material  to  be 
transferred.   Material  landed  on  the  transfer  station  of  the  combatant 
ship  is  cleared  from  the  area  by  fork  lift  or  pallet  truck.   Eventually 
the  transferred  material  will  be  transported  vertically  down  to  its 
designated  storeroom  or  magazine  by  elevator. 

The  type  of  movement  just  discussed  employed  the  unit  load  concept 
required  by  the  model  developed  herein.   As  noted  earlier,  many  small 
combatant  ships  must  break  up  the  unit  loads  into  smaller  units  for 
movement  beyond  the  transfer  station.   Many  times  hand  movement  is 
required  because  of  space  limitations  in  the  ship's  internal  passages. 
The  approximate  model  developed  by  this  study  could  still  be  used  for 


this  type  of  movement  using  a  mean  service  time  for  the  unit  load 
transferred  by  pieces. 

Queues  of  material  and  subsequent  blockages  of  material  can  and 
do  form  at  several  points  along  this  path.   On  the  supply  ship,  a 
queue  of  material  will  probably  exist  on  the  transfer  deck  near  the 
point  where  vertical  movement  from  the  hold  is  completed.   A  queue 
of  unit  loads  will  also  exist  at  the  transfer  station  as  material 
awaits  movement  to  the  combatant.   The  transfer  station  on  the 
receiving  ship  is  a  critical  blocking  point  as  transfer  of  a  load 
from  the  supply  ship  to  the  combatant  may  not  begin  until  the  landing 
area  has  been  cleared  of  the  previous  load.   Generally,  material 
cleared  from  the  landing  area  on  the  combatant  ship  is  temporarily 
stored  in  what  is  referred  to  as  a  marshalling  area.   Material  in 
the  marshalling  area  is  distributed  to  a  strikedown  route  and  then 
to  eventual  storage  in  the  assigned  storeroom.   In  many  cases, 
additional  queues  of  material  will  form  along  this  strikedown  route, 
especially  at  the  top  and  bottom  of  vertical  movement  leg(s). 

Much  nf  the  material  transferred  at  sea  presently  is  stored  as  an 
individual  component  or  packaged  with  similar  items  in  such  a  manner 
that  hand  movement  i*    p  These  individual  components  are 

combined  into  loads  that  can  be  conveniently  handled  by  the  transfer 
rig  between  the  two  ships.   This  load  makeup  is  accomplished  either 
in  the  hold  or  on  the  main  deck  of  the  supply  ship.   The  process  is 

entially  revet"     n  the  receiving  ship  with  the  breakup  of  loads 
commenced  early  in  the  storage  routine,  especially  on  smaller  ships, 
where  it  is  generally  cleared  from  the  landing  area  piece  by  piece. 
Introduction  of  more  and  better  handling  equipment  on  all  ships  and 
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movement  of  material  via  unit  loads  through  all  or  most  of  the  cycle 
will  vastly  improve  the  efficiency  of  the  operation.   This  step  forward, 
although  far  from  universal  reality,  is  in  progress. 

Careful  consideration  indicates  that  the  three  distinct  entities 
having  a  major  effect  on  the  movement  of  material  through  the  system 
are  the  service  or  handling  rates  through  all  phases,  the  storage 
capacity  of  the  various  queues  that  inherently  form  and  the  origination 
and  final  distribution  of  the  transferred  material.   The  importance  of 
the  handling  rates  is  obvious  to  the  extent  that  discussion  here  is  un- 
necessary,  However,  the  importance  of  the  queuing  (temporary  storage) 
space  and  the  distribution  of  material  requires  some  explanation.   Any 
amount  of  free  open  space  on  any  ship  is  always  scarce,  this  includes 
space  for  temporary  storage  of  transferred  material.   Perhaps  the  single 
exception  is  the  hanger  deck  on  an  aircraft  carrier,  (aircraft  carrier 
transfer  stations  are  located  off  the  hanger  deck),  however,  even  here 
space  may  be  at  a  premium  as  it  is  normally  taken  up  by  aircraft  in 
various  states  of  overhaul.   Therefore,  the  transfer  process  is  often 
slowed  because  of  blockages  caused  by  the  inherent  space  limitations 
on  board  ship. 

Limited  space  on  ships  is  also  partially  responsible  for  transfer 
slowdowns  that  can  be  attributed  to  excessive  distribution  to  a  single 
material  handling  facility,  e.g.,  a  vertical  strikedown  route  on  the 
combatant  ship.   Space  available  on  a  ship  restricts  both  the  size  of 
storerooms  and  vertical  and  horizontal  passages  leading  to  the  store- 
rooms.  In  order  to  provide  for  convenient  and  rapid  retrieval  of  a 
given  item,  it  must  be  stored  in  a  predesignated  place.   Further, 
many  items  require  special  storage  facilities  that  are  available  in 
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only  a  single  location.   These  special  facilities  can  range  from 
refrigerated  storage  for  fresh  and  frozen  provisions  to  storage  trays 
and  racks  for  most  ordnance.   Even  if  the  above  difficulties  could  be 
overcome,  space  limitations  again  would  be  limiting  in  that  storerooms 
are  generally  kept  full  so  that  the  only  storage  space  available  is 
the  designated  one.   Therefore,  it  is  often  impossible  to  change 
material  destination  to  accommodate  transfer  efficiency.   A  possible 
solution  to  blockage  on  a  route  would  be  an  alternate  route  for  some 
of  the  material.   Unfortunately,  alternate  routes  may  not  always  exist, 

B.   PARALLEL/SERIES  QUEUING  MODELS 

Solutions  for  a  single  stage  of  several  parallel  queues,  with 
various  types  of  input  and  service  times,  have  been  used  for  some  time 
Knowledge  of  single  server  queues  in  series  is  more  limited  and  is 
generally  restricted  to  Poisson  input  to  the  first  stage  and  exponen- 
tial service  times  in  all  stages  with  infinite  queues  allowed. 
Information  in  the  literature  of  a  combination  of  queues  in  series 
and  parallel  is  virtually  nonexistent. 

Saaty  [6J  gives  a  survey  in  the  area  of  series  queues.   The 
bibliography  for  this  work  is  nearly  complete  and  he  has  recently 
addended  the  bibliography  to  include  recent  contributions  in  the 
field   7  .   Burke   I   reported  the  fundamental  result  for  queues  in 
series  with  Poisson  input  to  the  fii      *go ,  exponential  service 
time  in  all  stages  with  infinite  queues  allowed,  showing  that  the 
output  from  each  stage  and  therefore  the  input  to  the  next  stage 
is  also  Poisson.   Hunt  [4J  using  the  assumption  of  exponential  service 
times  did  initial  work  in  the  area  of  allowing  only  finite  queues  to 
form  before  each  series  facility.   Hillier  and  Boling  [2J  extended 
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Hunt's  work  under  the  exponential  service  time  assumption  and  also 
solved  the  problem  for  Erlang  service  times  for  each  stage, 

Hillier  and  Boling  concurrently  reported  an  approximate  procedure 
for  determining  the  mean  output  rate  for  a  system  of  exponential 
service  time  facilities  in  series  with  finite  queues  allowed  before 
each  facility  except  the  initial  stage  where  an  infinite  queue  was 
assumed  to  exist.   Hillier  and  Boling  compared  results  obtained  with 
the  approximate  procedure  with  the  limited  exact  results  obtained. 
The  comparison  showed  very  little  error  in  all  cases  with  the  magnitude 
of  the  error  decreasing  with  increased  allowable  queue  size. 

The  Hillier  and  Boling  approximate  procedure  has  the  distinct 
advantage  of  being  computationally  feasible,  even  for  very  large 
problems,   This  approximate  procedure  formed  the  basis  for  the  study 
reported  here  and  it  will  be  discussed  and  outlined  in  some  detail  in 
the  next  section, 

C.   APPROXIMATE  PROCEDURE 

The  important  features  of  the  Hillier  and  Boling  f~2j  procedure, 
for  approximating  efficiently  the  mean  output  rate  (R)  for  a  series 
of  single  server  exponential  service  time  facilities,  may  be  outlined 
as  discussed  below.   Assume  that  there  are  N  facilities  in  series  each 
with  an  exponential  service  rate  of//,    i  =  1,2,  '"    ,  N,  and  waiting 
room  M.j_,  i  =  2,3,  '•■  ,  N,  where  the  M • ' s  include  the  item  being  served 
and  that  an  infinite  queue  always  exists  before  the  initial  facility. 

The  basis  for  the  procedure  comes  from  the  results  due  to  Burke  (_1J 
which  indicated  that  when  the  M^ ' s  are  equal  to  infinity  then  each 
service  facility  would  have  a  Poisson  input.   Reich  [5j  showed  further 
that  in  this  case  the  queue  sizes  would  be  independent.   This  implies 
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that  each  facility  could  be  analyzed  separately  using  the  standard 

M/M/l  model.   The  approximate  procedure  with  finite  queues  allowed 

is  then  based  on  the  assumption  that  the  output  from  each  stage  is 

approximately  Poisson.   Each  facility  is  analyzed  individually  by 

using  the  single-server  queuing  model  M/M/l  with  finite  queue  allowed. 

Note  that  R  is  merely  the  mean  output  rate  for  each  of  the  facilities. 

The  procedure  then  assumes  that  there  exists  an  effective  mean 

arrival  rate  A.(J'    (when  there  isn't  blocking  in  the  (j-1)    facility) 
ef  f 

to  the  j    facility.   It  is  also  assumed  that  the  queue  lengths  before 
each  facility  are  independent  so  that  the  unconditional  probability 
(1-pO-D)  that  the  (j-l)st  facility  is  busy  is  equal  to  the  conditional 
probability  that  it  is  busy,  given  that  the  queue  for  the  jtn  facility 
is  not  full.   Thus  it  follows  that, 

CD      *&\   =[^i-l(1-PoJ'U)  ■   «  n<M. 

^  °  ,   if  n  =  Mj  +  l 

where  n  is  the  number  of  customers  in  the  queue  for  the  jtn  facllll 

s  t 
(including  any  customer  being  held  at  the  (J-1)    facility  because  of 

blocking)  and  j  =  2,3,  ■••  .N. 

The  effective  service  rate  for  the  j   facility,  jtQffi   will  be 

somewhat  Lesfl  than  //•  .  j  =  2,3,  ••   .'  -1  ,  because  time  spent  in  holding 
a  customer  after  service  due  to  blockage  in  the  (j+l)st  queue  must  be 
included.   Now  the  mean  output  rate  of  a  single- server  queuing  system 
is  the  product  of  the  mean  service  rate  and  the  probability  that  the 
server  is  busy.   Therefore,  it  follows  that, 

(2)     y^ff  =  r/^-Pq-1^      •  for  j  =  2'3'  '"  'N- 

However,  for  the  final  server  the  customer  is  always  released  immedi- 
ately, i.e.,yY^l  =^"N.   Therefore, 
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(3)        R=^N(1-P^N)). 

Thus  R  may  be  determined  if  pW  is  known. 

o 

For  a  single  server  queuing  system  with  finite  queue,  it  can  be 
shown  that  [3 J 


y 


(1-y  MJ+1) 


(4)        }    p(j)  _  Jj H 

F°   "       M,+2 
l-yj  J 

where  for  the  case  of  the  approximate  system, 

It  is  also  clear  that  since  the  queue  before  the  first  facility  is 

V 

infinite,    P(1)  =  0  or  (1-P(1))  =  1. 
o  o 


It  is  now  possible  to  determine  the  true  value  of  R  by  first 

(2) 
*eff 


(1)  si  { 

assuming  a  trial  value  of  R,  call  it  R*.   Since  (1-PQ  )  is  known,  /\ 


can  be  determined  from  equation  (1) .   It  then  can  be  used  along  with  R* 
in  equations  (2)  and  (4)  to  determine  (1-P^2^).   Then  l-P^)  ,  (1-P^4^), 

•  •  •   ,  (1-P   )  can  be  obtained  in  a  similar  manner  so  that  R  can  be 

'      o 

calculated  from  (3).   This  calculated  value  of  R  is  compared  to  R*  and 
a  new  trial  value  assumed  and  the  process  repeated  until  the  true  value 
of  R  is  obtained  when  R*  -  R  =  0. 

Hillier  and  Boling  have  simplified  the  above  procedure  for  solution 
by  numerical  methods  and  the  equations  used  in  this  study  are  presented 
below: 

Combining  equations  (1)  and  (2)  yields, 

~#   (1-P<j"1))(l-P  (j))/R   ,   j  =  2,3,  •••  ,  N-l 
'   ij-1    o         o 

,   J  =  N 


yj  'K-^-eVA 
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then, 

(5)     i-p<^  =        R y,  ,  J  =  2,3,  •••  ,N-1, 

A-ia'*°    > 


(j-i)      J 


so  that  y.  can  now  be  found  by  combining  (4)  and  (5),  as  the  positive 
real  root  of  the  equation, 


M.  +  l 

(6)      d-y,J  )       n 

yj    J £ v  =  0 

K.+2         „   (J-D 


1"yjj      ^-l(1-Po    > 


Hillier  and  Boling  show  that  this  root  exists  and  is  unique  for  values 
of  R   that  are  sufficiently  close  to  R,  the  true  mean  output  rate. 
Let, 

M.  +  l 
) 

f  (v  )  =  "  J    -i     .  i  -  2'3»  '"  'N"1 

J   J     .     M.+2 


and 

g.(R,yj)  =  f.(y.)  -  c.y.,  j  =  2,3,  •••  ,N-1 
where  R  is  a  trial  value, 

-  _    R 


*    * 
and  X.  =  f.(y.),  >'.  corresponding  to  the  unique  root  of  g(R,y.)  =  0, 

ft 
Finally  X      vields    the   calculated   value   of   R   where 

y        N   N   ' 

S*  N 
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III.   FORMULATION  OF  THE  PROBLEM 

For  the  replenishment  at  sea  transfer  problem,  as  with  many  other 
physical  problems  that  can  be  formulated  in  a  queuing  theory  context, 
one  of  the  interesting  questions  that  can  be  asked  is,  "How  long  will 
the  operation  take?"  This  question  was  answered  in  the  formulation 
for  the  UNREP  transfer  problem  by  determining  the  mean  output  rate  (R) 
for  the  system. 

As  previously  noted,  each  at  sea  transfer  of  material  will  be 
somewhat  different  from  a  queuing  standpoint  because  of  the  ships 
involved,  volume  and  type  of  material  transferred,  material  origin 
and  destination,  types  of  material  handling  devices  employed,  etc. 
Thus,  the  number  and  size  of  queues  formed  during  any  transfer  may 
differ  considerably. 

A  completely  general  model  to  represent  the  transfer  at  sea  opera- 
tion would  require  up  to  five  parallel  paths  from  the  supply  ship  to 
the  combatant  ship.   On  each  of  the  parallel  paths  there  are  several 
material  handling  facilities  with  associated  queues  and  waiting  space 
on  both  the  sending  and  receiving  ships.   Additionally,  the  general 
model  would  require  movement  of  material  from  one  parallel  path  to 
another  at  least  once  on  each  ship. 

The  general  model  described  above  may  be  simplified  considerably 
and  still  represent  accurately  many  transfer  operations.   The  model 
developed  for  this  study  utilized  two  parallel  paths  from  the  hold 
of  the  supply  ship  to  a  marshalling  area  on  the  combatant  ship.   From 
each  of  the  two  marshalling  areas  the  material  is  moved  to  any  one  of 
three  strikedown  routes.   Figure  1  is  a  pictorial  representation  of 
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this  model.   The  inclusion  of  only  two  parallel  paths  was  done  for 
computational  efficiency,  noting  also  that  the  vast  majority  of  at  sea 
transfers  are  conducted  with  but  two  attached  rigs.   A  generalization 
of  the  model  for  more  than  two  paths  is  discussed  in  section  VI. 
Restriction  of  movement  between  parallel  paths  on  the  supply  ship  is 
the  other  major  simplification  made  for  this  study.   It  was  felt  that 
such  an  assumption  was  not  overly  restrictive  and  accurately  represents 
many  actual  transfers. 


Facility 
(11) 

11 


>* 


11 


Facil i 
(21) 
21 


Fac i 1  I 
(31) 
31 


M 


31 


Facility 

(12) 
12 


M 


21 


Facilitv 
(22) 
22 


M 


22 


Facil i  tv 
(32) 
32 

M 
32 


Facility 
(41) 
41 


M 


41 


Facil itv 
(42) 
42 


M 


42 


Fac  i 1 i tv 
(43) 
43 


M 


43 


Figure  1 

If  individual  material  handling  facilities  in  parallel  are  referred 
to  as  a  stage,  then  the  physical  location  of  the  first  stage  queue  of 
the  model  can  represent  the  hold  of  the  supply  ship.   The  second  stage 
queues  are  formed  by  material  on  the  main  deck  of  the  supply  ship 
awaiting  intraship  transfer.   The  marshalling  areas,  which  will  gener- 
ally be  adjacent  to  the  transfer  stations  on  the  receiving  ship,  form 
the  waiting  space  for  the  third  stage  queues.   The  fourth  itage  queues 
may  be  located  physically  at  any  one  of  several  places  on  the  individual 
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strikedown  routes  to  the  designated  storerooms,  e.g.,  on  the  transfer 
deck  near  an  elevator  or  hatch,  or  at  the  bottom  of  a  vertical  ^^vement 
leg. 

As  stated  previously,  service  of  material  in  the  two  third  stage 
queues  consists  of  movement  to  any  one  of  the  three  strikedown  routes. 
In  the  model,  this  movement  is  governed  by  two  probability  distribu- 
tions p^  and  q.,  i  =  1,2,3,  for  the  upper  and  lower  routes  respectively. 
The  distributions  are  determined  as  that  percent  of  the  material  moving 
across  one  of  the  initial  routes  that  is  scheduled  for  eventual  storage 
in  a  storeroom  on  the  i   strikedown  route.   Some  of  the  factors  that 
effect  this  determination  are  the  amount  and  type  of  material  trans- 
ferred, the  storage  of  material  on  the  supply  ship,  the  location  of 
the  storage  space  available  on  the  combatant  ship,  and  the  availability 
of  alternate  routes  to  final  storage.   It  is  felt  that  the  distribution 
approach  as  a  means  of  representing  the  actual  movement  during  transfer 
is  reasonably  flexible  to  cover  accurately  many  transfer  operations.   A 
primary  consideration  in  this  regard  is  the  fact  that  material  is 
selected  for  movement  in  the  most  convenient  manner,  i.e.,  a  load  is 
selected  for  movement  because  it  is  available  to  the  transfer  facility 
or  because  it  can  be  moved  ahead  of  other  unit  loads  to  a  nonblocked 
facility  upstream. 

Within  the  model  it  is  assumed  that  all  material  is  moved  as  a 
unit  load.   Solution  for  smaller  units  would  overly  complicate  the 
problem,  however,  the  unit  load  assumption  is  enhanced  by  the  present 
trend  toward  actual  movement  of  material  in  this  manner.   This  assump- 
tion had  a  further  and  probably  a  greater  impact  on  obtaining  an 
eventual  solution  because  of  the  effect  it  had  on  the  transfer  rates. 
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Since  all  material  moves  in  this  manner  it  is  reasonable  to  assume  that 
a  given  facility  can  handle  a  unit  load  of  500  lb.  bombs  or  1000  lb. 
bombs  or  even  fresh  provisions,  for  example,  at  approximately  the  same 
rate.   Thus  the  unit  load  concept  helps  justify  the  requirement  that 
each  facility  must  move  all  types  of  material  at  a  fixed  rate  during 
each  transfer  investigation.   There  is  of  course  some  allowance  for 
variation,  e.g.,  if  two  types  of  material  had  different  service  rates 
for  a  given  stage  and  each  moved  on  different  parallel  paths,  then  the 
service  rates  for  the  two  paths  could  be  adjusted  to  represent  more 
accurately  the  material  being  handled. 

Notation  within  the  model  is  such  that  the  first  subscript  (i) 
represents  the  i    stage  in  series  and  the  second  subscript  (j)  repre- 
sents the  j    parallel  element  in  a  stage,  e.g.,  y^*\2    *-a  tne  service 
rate  for  the  second  parallel  facility  in  the  first  stage  and  M^  is  the 
number  of  unit  loads  that  can  queue  before  the  first  parallel  facility 
in  the  third  stage. 

IV.   PARALLEL  SERIES  SOLUTION 

The  flow  of  material  through  the  physical  model  of  figure  1  may  be 
represented  as  in  figure  2.   If  R  is  the  mean  output  rate,  it  follows 


R 


II 


'III 


R 


Figure  2 
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from  previous  discussion  of  the  p.  and  q.  distribution  and 
figure  2  that, 

R  =  R.  +  R 

1     2 

R  -  RX  +  RXI  +  Rln 

RI  =  PlRl  +  «2*2 

R    =  p  R   +  q  R 
II    2  1     2  2 

RIII  "  P3R1  +  q3R2 

where  R  ,  R  ,  R  ,  RTT,  and  RJTT  are  the  mean  output  rates  of  their 
respective  routes.   For  any  UNREP  that  can  be  represented  by  this  model, 
the  value  of  the  mean  output  rates  will  be  determined  by  the  various 
service  rates,  queue  waiting  space,  and  the  p^  and  q^  distributions. 

If  the  system  were  symetric,  that  is  if  the  service  rates  and  queue 
lengths  for  the  facilities  in  each  stage  were  the  same  and  the  p.  and 
q.  distributions  moved  material  symetrically  to  the  final  stage,  then 
R  would  equal  R^ .   A  solution  to  this  problem  could  be  easily  reached 
by  the  Hillier  and  Boling  procedure  by  considering  only  half  the  system 
and  solving  for  R/2.   However,  these  conditions  will  not  often  hold  and 
in  general,  R  will  not  equal  R  . 

When  R  and  R  are  not  equal,  the  Hillier  and  Boling  approximate 
procedure  cannot  be  applied  directly  to  the  system  to  obtain  a  solution 
for  R.   The  problem  is  that  not  only  must  a  trial  value  for  R  be 
assumed,  but  also  values  for  R,  and  R  must  be  assumed  such  that 
R^  +  R~  =  R*,  where  the  R,r,s  are  trial  values.   Now  there  are  infi- 


nitely many  combinations  of  R  =  R*  +  R*  =  R,f.    The  problem  thus 


2  -  R  , 
reduces  to  finding  the  values  of  R^  and  Ro  that  will  yield  the  true 
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mean  output  rate.   The  true  values  of  the  R's  will  be  complicated 

functions  of  //     ,  M   ,  p  ,  and  q  .   It  will  not,  however,  be  necessar^ 
^Hj    ij    i       i 

to  determine  this  function,  as  the  solution  may  be  reached  by  simply 
maximizing  the  sum  R^  +  R2  for  a  given  value  of  R  . 

Before  proceeding  to  show  how  the  true  values  of  R, ,  R   and  R  may 
be  obtained,  the  Hillier  and  Boling  equations  should  be  adapted  to  the 
problem.   Assume  that  for  a  given  problem  the  true  values  of  R. ,  R  , 
and  R  are  known.   The  solution  for  the  output  across  the  upper  route 
may  be  found  by  applying  the  approximate  procedure  previously  outlined. 
First  >'   may  be  found  as  the  unique  positive  real  root  of 

g   (R  ,v   )  -  £   (v   )-c   v    =0 

621   1   21     21   21    2r  21 

where       c     R   /  <y  x"   ,  |  X      I).   Then, 

21     1  /   1 1  1  1      11 

(  1  -  v  1 1   ) 
r  21      21        ■ 

";1   =  21(y21}  =   i.v*M71+l 

l  ■ 

-k  if 

y    and  then  X   ma  h^     found  in  a  similar  manner.   Also,  X*  may  be 
found  using  the  iMM  technique  using  R*  in  place  of  R.  .   Then  since 

if 

the  fourth  stage  is  the  final  stage  X,.   i  ■  1,2,3,  can  be  determined 
direc  t 

X*   =  f   (     !X31Pi  *     'X32qi       i  =  1,2,3 
4i    4i  -JT 


so    that    finally   R   =   X*,^    ♦   X42^2      +  X' \jA 


ou 


43. 

The  solution  procedure  for  the  interesting  case  of  unequal  parallel 

v.-  & 

tputs  can  now  be  outlined.   As  before,  R   is  assumed  and  them  R^  and 


R2  are  chosen  such  that  R^  +  R2  =  R  .  The  problem  is  then  solved  as 
just  outlined  to  obtain  a  computed  value  of  the  mean  output  rate,  R. 
Using  the  same  value  of  R*,  R*  and  R2  are  varied  (maintaining  the 
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relationship  R£  +  Ro  =  R")  until  the  maximum  computed  value  of  R  is 

found.)   As  in  the  original  Hillier  and  Boling  procedure,  the  computed 

*  it 

value  of  R  is  compared  to  R'  and  if  R  =  R  ,  the  true  mean  output  rate 

has  been  found.   If  R  is  not  equal  to  R  ,  then  a  new  value  of  R   is 
assumed  and  the  procedure  repeated  until  the  true  mean  output  rate  is 
obtained.   This  procedure  converges  rapidly  if  two  techniques  are 
employed.   When  selecting  a  new  value  for  R'%  select  R  /new\=(Rx+R)/2. 
Also  when  commencing  each  iteration  R^ (new) /^2 (new)  =  ^T/^2. 

Justification  for  this  procedure  follows  from  physical  flow  laws. 
It  is  easy  to  see  that  the  procedure  is  correct  at  the  boundaries.   For 
example,  assume  that  the  system  was  such  that  one  of  the  //.  .'s  =  0  in 
the  lower  path.   Then  maximum  R  will  be  obtained  when  Rx  =  R-j"  and  R«  =  0. 

A  computer  program  to  solve  the  problem  of  figure  1  is  presented  on 
page  34.   The  program  was  employed  to  obtain  the  results  presented  in 
section  V.   The  wide  range  of  input  data  was  handled  with  no  difficulty. 
Average  computation  time  for  the  IBM  360  was  2  seconds. 

V .   PRESENTATION  OF  RESULTS 

As  noted  previously,  the  three  types  of  parameters  within  the 
approximate  transfer  model  that  have  an  effect  on  the  mean  output  rate 
are  the  service  rates  of  the  material  handling  facilities  (/^j^'s),  the 
available  waiting  room  in  each  queue  (M^'s),  and  the  distribution  of 
material  between  the  series  facilities  as  governed  by  the  p^  and  q^ 
distributions.   The  model  can  be  employed  to  study  UNREP  transfer 
operations  in  at  least  two  ways.   First  of  all,  general  parametric 
studies  can  be  made  to  determine  the  effect  of  varying  the  parameters 
ir  al '  possible  combinations.   Also,  more  specific  studies  may  be  made 
employing  known  handling  rates,  queue  spare  and  expected  material 
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distribution  for  specific  ship  types,  to  determine  critical  areas  that 
can  be  improved  to  increase  overall  transfer  efficiency.   A  second  use 
for  the  model  would  be  as  an  aid  to  the  commander,  perhaps  even  on  as 
low  a  level  as  commanding  officer  of  a  supply  ship,  to  help  him  in  the 
planning  for  replenishment  operations.   All  that  would  be  required  is 
a  knowledge  of  material  handling  characteristics  of  the  ships  involved, 
the  amount  of  material  to  be  transferred  and  a  small  computer  which  is 
generally  carried  on  board  many  supply  ships  today. 

The  primary  advantages  of  the  approximate  procedure  over  known 
computer  simulation  models  are  the  simplicity  of  the  analytical  model, 
the  inherent  computing  efficiency  and  the  ease  with  which  inputs  into 
a  computer  program  may  be  changed.   Possible  objections  to  the  model 
over  computer  simulation  models  would  include  the  limitation  that 
different  I      of  material  must  be  transferred  through  a  facility  at 
the  same  rate  and  that  unit  loads  of  each  type  of  material  would  occupy 
the  same  volume  or  area  of  queue  space.   The  exponential  service  time 
assumption  with  its  memorvless  prop*  | I     ...  the  probability  of  a 
unit  completing  service  in  a  time  interval  At  is  constant,  no  matter 
how  much  time  has  elapsod  since  the  last  completion  of  service,  also 
may  be  partially  objectionable.   A  final  objection  may  be  raised  over 
the  final  output  of  the  model  which  is  a  steady  state  result,  i.e.,  it 
is  assumed  that  the  system  has  been  in  operation  for  a  long  time  and 
is  in  equilibrium.   This  objection  is  at  least  partially  overcome  by 
the  generally  large  volume  of  material  transferred  during  the  UNREP's 
that  are  of  interest  and  the  comparatively  small  queue  sizes  allowed, 
thus  assuring  that  equilibrium  conditions  will  be  nearly  attained 
during  an  operation. 
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Actual  UNREP  transfer  data  was  not  available  in  the  format  required 
for  entry  in  the  model.   Therefore,  numerical  results  were  necessarily 
confined  to  several  parametric  studies.   Graphs  of  the  study  results, 
along  with  the  input  data  for  the  constant  parameters,  is  presented  in 
Appendix  A. 

Parametric  studies  1,  2  and  3  all  vary  the  parameter  p  ,  distribu- 
tion of  material  from  route  R1  to  route  RT ,  over  the  range  0.333  to  1.0. 

The  parametric  variation  in  the  three  studies  involves  the  q.  distribu- 
r  i 

tion  as  follows: 

Study  one  (page  29)  -  -  -  q.  =  1/3,   i  =  1,2,3 

Study  two  (page  30) qx  =  0.2,   q2  =  3=0-4 

Study  three  (page  31)  -  -  q  =  0.4,   q  =  q  =0.3 

All  three  studies  indicate  the  marked  decrease  in  overall  transfer 

rate  as  congestion  at  the  upper  strikedown  route  (facility  (4,1)) 

increases.   As  predictable,  the  decrease  is  most  rapid  for  study  3 

where  the  flow  to  facility  (4,1)  and,  therefore,  the  congestion  at  the 

facility  is  greatest.   An  additional  effect,  evident  in  all  three  studies, 

is  the  difference  between  R,  and  R   that  can  be  attributed  to  the  slight 

1      2 

reduction  in  handling  rates  of  facilities  (2,2)  and  (3,2)  and  the 
smaller  queue  allowed  at  facility  (3,2). 

Parametric  study  4  gives  a  good  indication  of  the  efficiency  that 
can  be  gained  by  increasing  the  handling  rates  of  the  facilities  in  a 
single  stage.   In  this  case,  the  service  rates  of  facilities  (2,1)  and 
(2,2)  were  varied  together  over  the  range  0.5  to  1.0  loads/min.  with 
all  other  service  rates  fixed  at  0.5  loads/min.   The  graph  of  the  study 
(page  32)  indicates  that  small  increases  above  the  system  average  will 
have  a  significant  effect.   Study  5  varied  the  queue  space  in  facilities 
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(3,1)  and  (3,2)  over  the  range  5  to  10  loads  with  the  remaining  queue 
spaces  held  fixed  at  5  unit  loads.   Once  again  a  graph  of  the  results 
(page  33)  indicates  that  small  changes  above  the  norm  are  worthwhile. 
Comparison  of  studies  4  and  5  shows  that  increased  service  rates  have 
a  greater  effect  on  efficiency  than  increases  in  allowable  queue  size. 
This  is  understandable  and  also  encouraging,  as  generally,  space  on 
ships  will  not  be  available  to  increase  temporary  storage  to  promote 
increased  transfer  efficiency. 

VI.   EXTENSION  OF  THE  MODEL 

There  are  several  wavs  that  the  specific  model  developed  here  can 
be  extended.   Unfortunately,  the  most  interesting  extension  becomes 
computationally  difficult. 

The  most  obvious  extension  follows  directly  from  the  Hi  1  Her  and 

Boling  \2~\    series  development  where  the  number  of  stages  in  series  is 

N.   The  parallel  s.ries  model  developed  here  may  be  similarly  extended. 

Referring  to  figure  2,  this  would  mean  that  there  could  be  N,  and  N 

1       2 

separate  service  facilities  in  the  two  initial  parallel  routes  and  N  , 
N   ,  and  N...  separate  service  facilities  in  the  three  final  routes, 
where  N   N  ,  N  ,  N   ,  and  N    are  only  restricted  to  be  finite. 
Solution  of  this  generalization  should  present  no  special  difficulty. 

The  interesting  extension  fron  an  UNREP  point  of  view  would  seem  to 
be  that  which  generalizes  the  number  of  parallel  paths  allowed.  Such  a 
generalization  is  represented  in  figure  3. 

In  solving  such  a  system,  it  is  first  of  all  evident  that  if  R  , 
R«,  •••  ,  R   are  known,  then  R,,  R,.,  •••  ,  RM  are  uniquelv  determined 
by  the  p      tri buttons .   Where  p..  is  the  distribution  of  material 
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from  the  i    initial  stage  to  the  j    final  stage.   Since  it  must  hold 
that 


and 


P±j  =  1,    i  -  1,2, 


j  =  1 
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Figure  3 
Also,  it  is  not  difficult  to  see  that  the  method  of  solution  presented 
for  n  =  2  earlier  becomes  increasingly  complex  as  n  increases.   The 
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true  computed  mean  output  rate  that  is  obtained  by  applying  the 

approximate  procedure,  is  found,  as  before,  by  determining  the 

-;.-  *        *    * 
combination  of  Ri+R~+  •••  +Rn  =  R   that  yields  maximum  R.   However, 

it  would  seem  that  an  iterative  procedure  could  be  developed  to  solve 

the  problem  efficiently  for  n  =  3  or  n  =  4  on  todays  high  speed 

computers.   Efficient  extension  beyond  this  point  would  probably 

require  knowledge  the  R.'s  explicitly  in  terms  of  the  service  rates, 

queue  lengths  and  flow  distributions. 
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M32  =  4  P2  =  P3  =  (1'pi)/2  qi  =  °'k  q2  =  q3 

PARAMETRIC  STUDY  3. 
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APPENDIX  A 
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~   .5  + 
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.*»-■ 


-  n 


O-  R  and  R2 

X-  RIt  Rn  and  Rni 


-o- 


-X- 


.5       .6      .7       .8      .9 
y^2\*  znd/Cf 22   (Second  Stage  Service  Rates) 


Values  of  other  parameters: 

Pl  -  qt  -  1/3,    i  «  1.2,3 


Xi-A2'A3  =  1/3 


PARAMETRIC  STUDY  *4. 
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APPENDIX  A 
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.8  4- 


.7  + 


.64 


.5 


.4.. 


.3- 


X 


.5 


-  R 


O  -  R-L  and  R2 

X-  Rx.  Rn  and  R^ 


-G- 


-*- 


-x- 


-x- 


.6 


.7 


.8 


-©- 


*7 


.9 


M-1  and  M_2  (Third  Stage  Queue  Length) 
Values  of  other  parameters: 

AmAwAz  =  1/3      VVV%  A3  =  5 

Pl  =  qi  =  1/3,   i  =  1,2,3 

PARAMETRIC  STUDY  5. 
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c 
c 


c 

c 


c 

c 

c 
c 
c 


c 
c 
c 
c 
c 


COMPUTER  PROGRAM  TO  OBTAIN  APPROXIMATE  SOLUTION  FOR  THE 
MEAN  OUTPUT  RATE  FOR  THE  REPLENISHMENT  AT  SEA  PROBLEM 

NOTE:  TO  SIMPLIFY  THE  COMPUTATIONS  THE  PROGRAM  HAS  BEEN 
WRITTEN  TO  SOLVE  ONLY  FOR  THE  CASE  WHERE  Rl  IS  GREATER 
THAN  OR  EQUAL  TO  R2.   THUS,  DATA  SHOULD  BE  ENTERED  TO 
MAKE  THIS  SO.   INSPECTION  OF  THE  HANDLING  RATES  FOR  REAL 
WORLD  DATA  WILL  NORMALLY  INDICATE  WHICH  PARALLEL  PATH 
WILL  HAVE  THE  LARGEST  MEAN  OUTPUT  RATE 

DIMENSIONED  VARIABLES 

YYU)  VALUES  ARE  USED  AS  THE  INITIAL  SOLUTION  FOR  THE 
RATIO  OF  THE  EFFECTIVE  ARRIVAL  RATE  TO  THE  EFFECTIVE 
SERVICE  RATE  VALUES.   VALUES  CORRESPOND  TO  1/10  OF  THE 
INTERVAL  (0,1)  AND  DEPEND  ON  THE  APPROXIMATE  MEAN 
SERVICE  RATE. 

M(I)  DEFINE  QUEUE  SPACE  AVAILABLE  AT  FACILITIES  (2,1) 
{2f2)t     (3,1),  (3,2),  (4,1),  (4,2),  (4,3). 

XS(I)  ARE  THE  PROBABILITIES  THAT  THE  FACILITIES  ARE 
BUSY  AND  IS  COMPUTED  WITHIN  THE  PROGRAM  FOR  FACILITIES 
(  l,l)-(2,3). 

U(I,J)  ARE  THF  EXPONENTIAL  SERVICE  TIME  INPUTS  FOR  THE 
FACILITIES.   NOTE  U ( 1 , 3 ) -U( 3, 3 )  ARE  DUMMIES  AND  ARE  NOT 
USED  IN  THE  COMPUTATIONS. 

P(I)  IS  THE  DISTRIBUTION  FROM  THE  UPPER  PATH  TO  THE  FINAL 
STAGE. 

Q(I)  IS  THE  DISTRIBUTION  FROM  THE  LOWER  PATH  TO  THE  FINAL 
STAGE. 

RHO(I)  IS  AN  INTERMEDIATE  VALUE  USED  IN  FINAL  COMPUTATION. 

RR(I)  ARE  THF  MEAN  OUTPUT  RATES  FROM  THE  FINAL  STAGE. 

DIMENSION  YY(10),M(7),XS(2,3),U(4,3),P(3),Q(3) 
DIMENSION  RHO(3),RR(3) 
WRITE(6,600) 
600  F0RMAT(T1 1  ,' WRITE  VALUE') 
RFAD(5,5)MM 
5  FCRMAT( 12 ) 
READ(4,46)YY 

46  FORMAT (  10F6.2) 
READ(4,47 )M 

47  F0RMAT(7I6) 
READ(4,48  )U 

48  F0RMAT(12F6.4) 
READ(4,49)P,Q 

49  FCRMAT(6F6.4) 
XS(  1,  1)  =  1. 
XS(  1,2)=1. 

1  =  1 
DEN=100. 

AR  IS  THE  INITIAL  TRIAL  SOLUTION  FOR  THE  MEAN  OUTPUT  RATE. 

AP=0.6*(U(4,  1)*U(4,2)«-U<4,3)) 
Rl  IS  THE  TRIAL  VALUE  FOR  THE  UPPER  ROUTE 

Rl=AR/2. 

R2  IS  THE  TRIAL  VALUE  FOR  THE  LOWER  ROUTE 
R2=AR-R1 

63  C=R1/U( 1,1) 
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SUBROUTINE  SELEC  DETERMINES  AN  INITIAL  TRIAL  SCLUTION  (Y) 
FROM  THE  TABLE  YY  FOR  THE  POLYNOMIAL  TO  BE  SOLVED  BY  THE 
NEWTON  RAPHSON  METHCO. 

CALL  SELEC(C,YY, Y) 


C 
C 
C 


SUBROUTINE  APROX  SOLVES  THE  POLYNOMIAL  TO  DETERMINE  THE 
RATIO  OF  EFFECTIVE  ARRIVAL  RATE  TO  EFFECTIVE  SERVICE  RATE 
(Y)  AND  THEN  COMPUTES  XS(I,J). 


60 


CALL 
C  =  R2/ 
CALL 
CALL 
C  =  R1/ 
CALL 
CALL 
C  =  R2/ 
CALL 
CALL 
A  =  XS( 
B  =  XS( 
DC    60 
RHOIK 
D=RHO 
RR(K) 


M( 1) ,Y,XSI2,  1)  ) 


2,2)  ) 


1)  ) 


APROXIC 
U(l,2) 

SELEC<C,YY,Y) 
APR0X(C,M(2) ,Y.XS( 
(U(2,1)*XS<2,1)  ) 
SELEC(C,YY,Y) 
APR0X(C,M(3) 
(U(2,2)*XS(2 
SELEC(C,YY,Y) 
APROX  (C,MH)  ,Y,X 
3fU*U(3,l  ) 
3,2)*U( 3,2) 

K=  1     3 
)=<P<  K)**  +  UK»*B)/U(4,K) 
(K)**(M(4+K) +1) 
=(RHO(K)-D)/ (l.-D)nU4,K) 


Y,XS(3, 
2)  ) 


c    ( 


UZ)  ) 


RS     IS    THE    COMPUTED    VALUE    OF     THE    OUTPUT    RATE. 

RS=RR< 1)+RR12)+RR(3) 

THIS  SECTION  OF  THE  PROGRAM  COMPARES  THE  TRIAL  VALUE  AR 
TO  THE  COMPUTED  VALUE  RS  AND  ITERATES  UNTIL  THE  DIFFERENCE 
IS  LESS  THAN  .001.   THE  SECTION  ALSO  CAUSES  ITERATIONS  TO 
BE  MADE  TO  FIND  THE  VALUES  CF  Rl  AND  R2  THAT  YIELD  MAX  RS . 


64 


61 


62 


65 


69 


70 


67 

68 


IF 
IF 

RS 
Rl 
R2 
GC 
RS 
1  = 
GO 
Rl 
R2 
IF 
DE 
1  = 
GC 
1  = 
DE 
S  = 
T  = 
IF 
IF 
AR 
IF 
R3 
IF 
Rl 
R2 
GC 
AR 
GO 
Rl 


I.EQ.l  )C()  TO  6  1 

RS.LT.RS1 )G0  TO  62 

=  RS 

Rl+l./DFN 

AR-R1 

TO  63 

=  RS 

TO  64 

Rl-l./DEN 

AR-R1 

I  .NE.2)G0  TO  65 

=1000. 


TO  64 

=  100. 

1/R2 

BS(RS-AR) 

T.LE.0.002)G0  TO  66 

RS.GT.ARJGO  TO  67 

AR-T/2. 

R1.E0.R2)G0  TO  68 

ABSIR1-R2) 

R3-.l)68,68,70 

AR/2.+R3-.1 

AR-R1 

TO  63 

AR+T/2. 

TO  69 

AR/2. 
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66 
55 


52 
53 
54 


R2=AP-R1 
GC  TC  63 


UU  IU  63 

WRITE(MM,55)RS,R1 ,R2,RR 

FCRMATI//T5, 'QUEUE  OUTPUT  IS',T23,F8. 

OUTPUTS  ARE',T28,2F8.5//T5, 'TERMINAL 

T28,3F8.5) 


5//T5, 'PARALLEL 
OUTPUTS  ARE', 


SLBROUTINE  APROX ( V AL , N , Y , X > 

AN  =  N 
30  A=Y**N 

B  =  A*Y 

C=( Y-B  )/(  l.-B)-VAL*Y 

D=(1.-(AN*1.)*A+AN*B)/(1.-B)**2-VAL 

Y1=Y-<C/D) 

IF(ABS(Y1-Y)-.000  1)20,20,22 
22    Y=Y1 

GC  TO  3C 
20  X=(Y-B)/< l.-fl) 

RETURN 

END 


25 
23 


C 

c 


SL 
DI 
DO 
IF 
CC 
B  = 
RF 


BROUTINE 
MFNS  ION 

25  1=1, 
(C-FLOAT 
NTINUE 
A<  I  ) 
TURN 
D 


SELECl 
A(  10) 

n 

(  I  )/  10. 


C,A,B  ) 


)23, 23,25 


SAMPLE  DATA  DECK 


YY(  I  ) 

NOMIAL 
YY(  1  )  YY 

20.01 
C 


M(  I 


I 


) 

M21 

5 


INITIAL  TRIAL  SOLUTION  FOR  NEWTON  RAPHSCN  POLY- 
SOLUTION,    IT  OEPENDS  ON  THE  VALUE  OF  AR  ASSUMED. 
(2)  YY(3)  YY(4)  YYI5)  YY(6)  YY(7)  YY(8)  YY ( 9 )  YY(iO) 
6.00  <*.00       2.80   2.20   1.80   1.5}   1.20   0.90   0.60 

MAX  QUEUE  LENGTHS 

M22     '  II     M32     M4  1     M42    M43 
5       5       4       5       5       5 


I 


U( 

ON 

Ul 

5000 

U33 

5333 

P(  I  ) 
P(  1) 
3333 


,J) 
A  S 


-  HANDLING  RATFS  (NOTE  FUR  ACTUAL 
INGLE  CARD  WITH  12F6.3  FORMAT). 

U31    U4l    U12    U22    U31 
5000  .3333  .5000  .4000  .4000 


INPUT  ALL  DATA  IS 


U2  1 
5000 
UA3 
3333 


U4  2 
3333 


U13 
50C0 


U23 

5000 


iND  Q( I )  DISTRIBUTIONS 
P(2)   P( 3)   Q( I )   Q(2)   Q( 3) 
3333  .3333  .3333  .3333  .3333 
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